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Abstract 

Light interacts differently with left and right handed three dimensional chiral objects, like helices, 
and this leads to the phenomenon known as optical activity. Here, by applying a polarization 
tomography, we show experimentally, for the first time in the visible domain, that chirality has a 
different optical manifestation for twisted planar nanostructured metallic objects acting as isolated 
chiral metaobjects. Our analysis demonstrate how surface plasmons, which are lossy bidimensional 
electromagnetic waves propagating on top of the structure, can delocalize light information in the 
just precise way for giving rise to this subtle effect. 
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I. INTRODUCTION 



Since the historical work of Arago [l| and Pasteur (2|, chirahty (the handedness of 
nature) has generally been associated with optical activity, that is the rotation of the plane 
of polarisation of light passing through a medium lacking mirror symmetry 0, 4|. Optical 
activity is nowadays a very powerful probes of structural chirality in varieties of system. 
However, two-dimensional chiral structures, such as planar molecules, were not expected 
to display any chiral characteristics since simply turning the object around leads to the 
opposite handedness (we remind that a planar structure is chiral if it can not be brought 
into congruence with its mirror image unless it is lifted from the plane). This fundamental 
notion was recently challenged in a pioneering study where it was shown that chirality has 
a distinct signature from optical activity when electromagnetic waves interact with a 2D 
chiral structure and that the handedness can be recognized {sj. While the experimental 
demonstration was achieved in the giga-Hertz (mm) range for extended 2D structures, the 
question remained whether this could be achieved in the optical range since the laws of 
optics are not simply scalable when downsizing to the nanometer level. Here we report 
genuine optical planar chirality for a single subwavelength hole surrounded by left and right 
handed Archimedian spirals milled in a metallic film. Key to this finding is the involvement 
of surface plasmons, lossy electromagnetic waves at the metal surfaces, and the associated 
planar spatial dispersion . Our results reveal how, in a stringent and unusual way, this 
optical phenomenon connects concepts of chirality, reciprocity and broken time symmetry. 

We remind that partly boosted by practical motivations, such as the quest of negative 
refractive lenses js] or the possibility to obtain giant optical activity for applications in 
optoelectronics, there is currently a renewed interest js- 18| in the optical activity in artificial 
photonic media with planar chiral structures. It was shown for instance that planar gam- 
madionic structures, which have by definition no axis of reflection but a four-fold rotational 
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can generate optical activity with giant gyrotropic factors [12 



16 



-ll8|. 



Importantly, and in contrast to the usual three dimensional (3D) chiral medium (like quartz 



and its helicoidal structure 



19| ) , planar chiral structures chaiige their observed handedness 



when the direction of light is reversed through the system [9|, |20[. This challenged Lorentz 
principle of reciprocity (which is known to hold for any linear non magneto-optical media) 



and stirred up considerable debate 



12 



2l| which came to the conclusion that optical 
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activity cannot be a purely 2D effect and always requires a small dissymmetry between the 



two sides of the system 
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16l4l8j. Nevertheless Zheludev and colleagues did demonstrate 



in the GHz spectrum that a pure 2D chiral structure lacking rotational symmetry can have 
an optical signature which is distinct from optical activity [5]. They went on to predict 
that it should be possible to observe the same phenomena in the optical range by scaling 
down their fish-scale structure and playing on localized plasmons 22|. Following a different 
strategy, we show here that SP waves propagating on a 2D metal chiral grating resonantly 
excited by light provide an elegant solution to generate planar optical chirality in the visible. 



II. EXPERIMENTS AND RESULTS 

This is a challenging issue as it leads to two fundamental points which are apparently 
incompatible. On the one hand, finding such a 2D chiral effect in the optical domain is 
not equivalent to a simple rescaling of the problem from the GHz to the visible part of 
the spectrum. Indeed, losses in metal become predominant at the nanometer scale so that 
the penetration length of light through any chiral structure will become comparable to the 
thickness of the structure. In-depth spatial dispersion along the propagation direction of 



light will hence be induced, corresponding to the usual 3D optical activity |lll-ll8|. One thus 
expects optical activity, through the losses, to be a more favorable channel than 2D optical 
chirality. On the other hand, losses (i.e., broken time invariance at the macroscopic scale) 
are necessary to guarantee planar chiral behavior j^, |22|. With this in mind, we chose to 
make single SP structures such as a single hole in an optically thick metal film surrounded 
by an Archimedian spirals (Figure 1) which can provide all the necessary ingredients for ob- 
serving 2D optical chirality. It is a 2D structure lacking point symmetry, that is rotational 
and mirror invariances. At the same time, it resonates due to coupling to surface plasmons 
which, as lossy waves, represent a natural way for delocalizing information along a planar 
interface, moving in-depth losses to the surface. Importantly, the thickness of the metal 
film optically decouples both interface ^|, and consequently only the structured chiral side 
is involved in the 2D optical chiral effect reported here. Finally, the structures gives rise 
to enhanced transmission 7| enabling high optical throughput for all the characterization 
measurements. 
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Using focus ion beam (FIB), we milled in an opaque gold film a clockwise (right TZ) or 
anticlockwise (left C) Archimedian spiral grooves around a central subwavelength hole. The 
polar equation {p,0) of the left handed Archimedian spiral is p = P ■ 6/{2tt), and the right 
handed enantiomeric spiral is obtained by reflection across the y axis (see Fig. 1). The 




Wavelength (nm) 

FIG. 1: Chiral plasmonic metamolecules. On the top panel: scanning electron micrographs of 
the left (L) and right (R) handed enantiomer (mirror image) planar chiral structures investigated. 
The scale bare is 3 i-im long. The parameters characterizing the structure are the following: hole 
diameter d = 350 nm, film thickness h = 310 nm, grating period P = 760 nm, groove width 
w = 370 nm, and groove depth s = 80 nm. The structures are milled, with a focus ion beam, in 
a gold film deposited on a glass substrate. On the bottom panel: transmission spectra at normal 
incidence of individual left (blue curve) and right handed (red curve) Archimede spirals illuminated 
from the air side. 

geometrical parameter P is the radial grating period and we take its value equal to the 
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SP wavelength Xspp — 760 nm (for an excitation at A ~ 780 nm). We recorded optical 
transmission spectra at normal incidence with unpolarized light for both isolated structures 
(Fig. 1). As it can be seen, both enantiomers behave like resonant antennas with quasi 
identical transmission properties. This resonant behaviour is a direct indication of the SP 
excitation by the grating similarly to what is observed for circular antennas (24] . 

To observe and fully characterize the optical signature of planar chirality we perform a 



full polarization tomography 



25 



26| with the aim of determining the 4x4 Mueller matrix 
Ai associated with each enantiomer. Experimental results A^^^'^P', and Ai-R,^'^'^' respec- 
tively obtained for left and right handed spirals are given in appendixes A and B. Here the 
important point is that the degree of purity F of the Mueller matrices 25|| is near unity with 
F (^A^£*^^P-) ~ 0.967 and F (A^7^^^P■) ~ 0.939. This shows that the coherence in polar- 
ization is not degraded by the structure and that we can therefore restrict our discussion to 



Jones matrices 



25| . In the convenient left \L) and right |-R) circular polarization basis, 
these Jones matrices tie the excitation [i?}^'^, to the transmitted [£'2^^, E"^^^] electric 
fields. In the case of planar chiral structures displaying 2D chiral activity, they have the 



following form [5|, l22| : 
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where A, B and C are complex valued numbers such that \B\ ^ \C\. This inequality 
account for chirality. Being non diagonal, these matrices correspond to polarization con- 
verter elements with no rotational invariance around the z axis (Fig. 1). They are thus 
fundamentally different from Jones matrices associated with optical activity, e.g., gamma- 
dions. Importantly the conditions \B\ ^ \C\ implies the non unitarity of JT^^ which means 
that reversing the light path through the chiral structures is not equivalent to reversing the 
time. From equation (1) we deduce the associated theoretical forms for the Mueller matrices 

f h tVi 

■J^c Ai-ji ' (see appendix C) which are used to fit Jc ^"^^ <Jn from experimental results. 



After normalization by A we deduce 

fit ' 1-000 



Jc 



Jn 



fit 



0.166 + i0.221 
-0.131 + i0.099 1.000 

1.000 -0.129 + i0.098 
0.170 + i0.230 1.000 



(2) 



These matrices indeed satisfy the chirahty criteria of equation (1) within the ~ 1% uncer- 
tainty evaluated from the degree of purity of the Mueller matrice of the empty setup. 




angle (rad.) angle (rad.) 

FIG. 2: Analysis of the polarization states for an input light with variable linear polarization for 
both the left (left panel) and right handed (right panel) individual chiral structures of Fig. 1. 
The data points (acquired with a laser light at A = 780 nm) are compared to the predictions 
from equation (2) (continuous curves) for respectively the transmitted intensity analyzed along 
the direction: |a;) (green), \y) (yellow), | + 45°) (cyan), | — 45°) (magenta), |L) (red), and \R) 
(blue). The total transmitted intensity is also shown (black). The symmetries between both panel 
expected from group theory (see appendix D) are observed experimentally. The insets show in each 
panel the ellipses of polarization and the handedness (arrow) associated with the two corotating 
eingenstates associated with the Jones matrix Jc (blue) and J-ji (red). 



III. DISCUSSION AND CONCLUSION 

To illustrate the polarization conversion properties of our chiral structures, we compare in 
Fig. 2 theory and experiment when the input state is linearly polarized and when the output 
transmitted intensity is analyzed along different orthogonal directions. A good agreement 
between the measurements and the theoretical predictions deduced from the Jones matri- 
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ces (see appendix D) is clearly seen, together with the mirror symmetries between the two 
enantiomers. Importantly, these symmetries also imply that for unpolarized light, and in 
complete consistency with Fig. 1, the total intensity transmitted by the structures is in- 
dependent of the chosen enantiomer. Furthermore, the conversion of polarization is well 
(geometrically) illustrated by using the Poincare sphere representation 25|. Indeed, as 
shown in Fig. 3, the Mueller matrix defines a geometrical transformation which projects the 
unit Poincare sphere, drawn by the input Stokes vector, on an output closed surface with 
typical radius F (A^^^P') ^ 1 in agreement with the absence of net depolarization as already 
noticed (from theory, F (^Ai'"^'^ = 1 exactly). Data shown on Fig. 2 are also plotted on 
this sphere. The input state draws a circle in the equator plane while the output state (for 
each enantiomer) draws a circle in a different plane, which center is not located at the center 
of the sphere. This is a direct manifestation of planar chirality (see appendix E). There 
is clearly an antisymmetrical behaviour between both enantiomers. The good agreement 
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FIG. 3: Full polarization tomography. Poincare sphere of unit radius associated with the input 
state represented by the Stokes vector X j^, [2^. Also shown are the results of Fig. 2 for the left 
(blue) and right handed (red) structures if the linearly polarized incident state draw the black circle 
in the {Xi, X2) equator plane of the input sphere. Data points are compared with the predictions 
from AicTl^^"^' (continuous curves) and of equation (2) (dashed curves). 

between the experiment and the prediction of equations (1,2) shows the sensitivity of the 
polarization tomography method and the high reliability of the FIB fabrication. 
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th th 

The degree of optical 2D chirality is quantified by diagonalizing Jc. ' and Jti '. 

-J-l-. y , 

For J'c ■, the eigenstates are \±c) = \B\L) ± y C\R) associated with the eigenvalues 
A£(±) = y4 ± ^{B ■ C). The eigenstates for Jti"^' are obtained by permutation of B and 
C with consequently A£(±) = A7^(±). The scalar product {+c\—c)^^' = —{+ii\—n)^^' = 
— |C|)/(|i?| + |C|) is the eigenstates Stokes parameter S3/S0 and provides a direct mea- 
surement of the degree of optical chirality. It also evaluates losses since the non-orthogonality 
of these two states is related to the necessary non-unitarity of the Jones matrix for planar 
chirality. We have {+c\—c)^^ — 0.255 and {+n\—n)^^ — —0.277, which, within experi- 
mental uncertainties, are in good agrement with the theoretical expectations. As shown in 
the insets of Fig. 2, both eigenstates of each structure (e.g., |±£,7^)) can be represented by 
two ellipses having the same axis ratio and the same handedness, but rotated 90° relative to 
each other. In agreement with the theoretical predictions, these polarization ellipses for |±£) 
and |±7e) are mirror reflections. This behaviour is significantly different from the results 



obtained with optically active media 



12-[l4l |l6;18| where the eigenstates associated with 



a given enantiomer have opposite handedness 4- This point, which reflects itself in the 
symmetry property of chiral Jones matrices, namely J'^\{L, L) = J^cni^^ ^) ~ ^1 
reaching consequences, as pointed out in reference [Sj . It implies that a 2D plasmonic spiral 
mimics a Faraday medium when we reverse the light path and this even if the system, unlike 
a true Faraday medium, obeys rigorously to the principle of reciprocity ji, 5| (inversely, one 
can show that equation (1) results from both this requirement and the absence of mirror 
symmetry). It means that a photon coming from the second side will probe a structure of 
opposite chirality. After going through the structure and retracing back the light path with 
a mirror normal to the axis, the polarization state will be different at the end of journey 
from the initial one. This would be impossible for an optically active medium and is solely 
due to planar chirality. To summarize, our results therefore demonstrate that 2D chirality 
is possible in the visible domain in the absence of optical activity and add another element 
to the promising plasmonic toolkit. 



IV. APPENDIX A: POLARIZATION TOMOGRAPHY SETUP. 



We apply a procedure similar to the one considered in [26|, |27[ in order to record the 
Mueller matrix: a collimated laser beam at A = 785 nm is focussed normally on the struc- 
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ture by using an objective Li (x50, numerical aperture=0.55). The transmitted light is 
collected and recoUimated by using a second objective L2 (x40, numerical aperture=0.6). 
The input and output states of polarization are respectively prepared and analyzed in the 
coUimated part of the light path by using polarizers, half wave plates and quarter wave- 
plates. A sketch of the setup is provided below (see Fig. 4). 
The Mueller matrix is built by applying an experimental algorithm equivalent to the one 




FIG. 4: Principle of the polarization experiment, (a), Sketch of the optical set up described in the 
text. The images are recorded by using a CMOS camera, (b), A typical image of the transmitting 
nanohole showing the Airy spot associated with diffraction by the optical microscope. The scale 
bar is 2 /im long, (c), Crosscut of the intensity profile along the yellow dotted line shown in (b). 



described in 25|. More precisely, in order to write down the full Mueller matrix, we mea- 
sured here 6x6 intensity projections corresponding to the 6 unit vectors \x), \y), | + 45°), 
I —45°), |L), and |-R) for the input and the output polarizations. Actually only 16 measures 
are needed to determine M. [25| . Our actual procedure is thus more than sufficient to obtain 
M. 

The isotropy of the setup was first checked by measuring the Mueller matrix Jv[^^^^ with 
a glass substrate. Up to a normalization constant, we deduced that A^S^^^^ is practically 
identical to the identity matrix X with individuals elements deviating by no more than 0.02. 
More precisely, the optical depolarization (i. e, the losses in polarization coherence) can 



be precisely quantified through the degree of purity of the Mueller matrix defined by 



25| 
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1 /2 

F (M) = {^ ^"^^^smI'^"" ) ^ 1- Here we measured F (^M^^^^^) = 0.9851. It implies that 
the hght is not depolarized when going through the setup and that consequently we can rely 
on our measurement procedure for obtaining Ai. 

Two important points must be noted here: On the one hand we varied the incident illumina- 
tion spot size on the sample between 2 and 20 fim without affecting the matrix, i. e., without 
introducing additional depolarisation. In the rest of the experiment on chiral structures we 
consider the case of a large gaussian spot with FWHM=20 yum in order to illuminate the 
whole individual spiral. On the other hand, it can be observed that in our experiments the 
polarization in the Airy spot (see Fig. 4b) is homogeneous. This implies that we are actu- 
ally doing the polarization tomography of the central transmitting hole, i. e., we are dealing 
only with the SU(2) point symmetry of the Mueller Matrix. This situation clearly contrasts 
with previous SOP tomography measurements on metallic hole arrays in which the polar- 
ization degrees of freedom were mixed with spatial information responsible for SPP-induced 
depolarization |26 |. 



V. APPENDIX B: EXPERIMENTAL MUELLER MATRICES. 



The experimental Mueller matrices deduced from the polarization tomography are after 
normalization of A^qo'^ ■ 



>^.exp. 



1.000 


0.031 - 


0.029 


0.958 


-0.105 


0.037 


0.029 


0.261 - 


1.000 


0.035 


0.027 


0.949 


0.096 


-0.034 


-0.011 


-0.252 



0.107 -0.029 



0.287 



0.111 0.023 



-0.277 0.745 



(3) 



We have F (tW^^^P') ~ 0.967 and F (A^Te^^P') ~ 0.939. 

We must also note that the normalization used here neglects a small additional coefficient 
of proportionality {Aicff^^' / -M-Tiff^^' \ — 0.954 imputed to experimental errors and uncer- 
tainties. 
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We also recorded the Mueller matrix of the set up with the glass substrate only. Up to a 
normalization factor we deduced 



1.0000 0.0060 
-0.0030 0.9851 



-0.0040 -0.0070 
-0.0010 0.0020 



-0.0020 0.0020 0.9965 0.0030 
-0.0050 -0.0040 0.0030 0.9821 



(4) 



which satisfies A4§^^^^ ~ X with X the identity matrix. It imphes that the optical set up do 
not induce depolarization and that consequently we can rely on our measurement procedure 
for obtaining M.. 



VI. APPENDIX C: THEORETICAL MUELLER MATRICES. 



th 

The precise form of the theoretical Mueller matrice M.c ' deduced from equation (1) is 
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(5) 



with M^^- = (2|Ap + \B\^ + |Cn/2, TWjf • = Re[BA* + AC*], XjJ- = Im[AB* + CA*], 
A^th. ^ Q^|2 _ |5|2)/2, m\^- = l^p + Re[B*C], M\^- = Im[B*C], M\^- = Re[AC* - 

BA*], mI^- = l^p - Re[B*C], M\}- = Re[CA* - AB*], M^}- = (2|A|2 - |Sp - |Cp)/2. 

th 

Similar formula are obtained for M.-R, ' after permuting B and C. 

th 

Prom the previous relations we deduce the useful equations (valid for A^^ ') 



B/A-- 
C/A 



A^th. 

•'^102 



th. 

13 



(6) 



Together with equation (3) equation (6) allow us to fit B/A and C/A if we replace JH^^- 
by A^£^^P" (a similar procedure is apphcable to M-r,^^^' after permuting B and C). 
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The best fit we obtained (see equation (2)) are: 

/ 



M 



fit 



7^ 



1.000 0.033 -0.116 -0.023 

0.033 0.951 0.043 -0.282 

-0.116 0.043 0.951 0.304 

^ 0.023 0.282 -0.304 0.902 

^ 1.000 0.0359 0.125 0.026^ 

0.039 0.949 -0.044 0.283 

0.125 -0.044 0.948 0.311 

-0.026 -0.283 -0.311 0.897 j 



(7) 



From theory we can deduce that F yM.c,n^j — 1 (i-e., after normalization by M^')- We 
have thus F (^Mcn^^^ = 1 



VII. APPENDIX D: SYMMETRIES DUE TO CHIRALITY [INTERPRETING 
FIGURE 2]. 

Let ~ Ex\x) + Ey\y) and |^out) ~ ^'x\^) + E'y\y) be respectively the incident and 

transmitted electric fields when we consider the left handed planar chiral structure. We have 



l^out) = M"^ 



ml 



(8) 



where Jc is the operator associated with the Jones matrix Jc- The mathematical definition 
of planar chirality is that whatever the mirror symmetry operation IT in the plane X-Y we 
have jfn — lij ^ 0. It equivalently states that tij'flr'^ ^ J . If we consider for example 
the mirror refiection through the Y axis (see Fig. 1) we have the matrix representation (in 

-1 
1 



the cartesian basis) 11 = 11 ^ 



and consequently 



(9) 



which agrees with equation (1) and constitutes an other optical definition of chirality. 
The previous equations are used in order to interpret the results of Fig. 3 of the main article. 
Indeed from equations (8) and (9) we obtain 



n|*out) = ^7en|*in). 



(10) 
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The input state considered in Fig. 2 is a linearly polarized light |^) = sin {0)\x) + cos {9)\y) 
(the angle is measured relatively to the Y axis) and the transmitted intensity projected 
along a direction of analysis \i) (i.e, \x), \y), \ + 45°), | — 45°), \L), and \R)) is written 
^(Left)(^) = K^l^out)!' = mJclO)?. Similarly we also write ll^'^^^\9) = mJn\9)?. 
From equation (10) we deduce: 



{^'\Jc\0) = {^\Jn\-0), (11) 
where we used = fl~^\i) = fl\i) and \ — 0) = fl\0). We consequently have: 

.(Left)..^ _ .(Right), n^ 

-'±45° — -'+45° \~^)-: 

Such symmetries are clearly visible in Fig. 2 and correspond to a direct signature of optical 
chirality in the planar systems considered. 

VIII. APPENDIX E: PLANAR CHIRALITY ON THE POINCARE SPHERE [IN- 
TERPRETING FIGURE 3] 

We remind that the Stokes parameters associated with a polarization state of light |\1/) 
are defined by 

So = Ix ~^ lyi Si = Ix — ly 
S2 = 1+45° - I~45°, 'S's = /l - /j^, (13) 

where Jj are projection measurement along the direction i, i.e, = |(z|^)p. The Stokes 
vector X is a convenient representation of such a state. We have X = XiXi + X2X2 + X3X3 
with Xi = Si/ So, X2 = S2/S0, X3 = S3/S0 and with (xi, X2, X3) a cartesian orthogonal 
and normalized vector basis. 



The coherent input state satisfies the normalization 251] |X| = 1, that is the vector draw a 



Poincare sphere of unit radius in the space Xi, X2, X^. The transmitted output state after 
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interaction with the left or right handed structure is defined by the relation 







( ^ \ 

Do 




= Mc,n 


Si 






S2 


\ Sc,n;3 ) 







(14) 



The output state defines a Stokes vector X^^t^ such that |X£^7^| < 1. A typical value for this 
radius is given by F {A4c,Ti)- 

If the input state is linearly polarized the input Stokes vector is: 

^ cos (2^) \ 
sin (2^) 




V 



(15) 



and draw a circle (Z^jn) along the equator contained in the plane Xi, X2 of the unit radius 
Poincare sphere. Using equation (14) the output Stokes vector is now a function of 9: 
^c,n{9) drawing a closed curve {Y^c^-ji) (see Fig. 3) which is the image, through the Mueller 
matrix transformation, of the equator circle (Z^j^) above mentioned. Importantly, since 
the Mueller matrix Ai given by equation (5) represents a linear relation connecting X^^ to 
Ij^out^ we conclude that the image of the incident polarization state contained in the equator 
plane Xi, X2 through A4 must also be contained in a plane in the space Xi, X2, X^. 
To analyze this point more in details we consider the normalized Vector product 

(X£,7^(0) - X^,^(27r/3)) x (X£,^(0) - X£,^(7r/2)) 



{Xc,n{0) - x (X^,7^(0) - XcM^/2))\ 



(16) 



and we write it 



^ Uc,n ^ 



^c,n — 



(17) 



with \Uc,Ti\'^ + |Vc,7^P + iVFc^T^p — 1. It represents a typical normal to the closed curve 
(E£,7^)- We have 



0.2845 
-0.3065 
-0.9084 



0.2861 
0.3139 
0.95053 



(18) 
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Actually, if each curve {J2c,tz} is contained in a (different) plane Pc,iz we must have 

nc,n ■ (Xr,7^(^) - X^,7^(0)) = (19) 

for every 9. This was indeed checked numerically up to a precision of 10~^^. It was also 
checked that |X£^7^(^)| = 1 up to the same precision. This proves that each curve (J^cr) 
must be a circle. The equations of the two planes Pc,n are given by n^ T^- (X — X£ 7^(0)) = 
where X is the Stokes vector associated with a running point belonging to each plane. We 
write 

Uc,nX, + Vc,nX2 + Wc,nXs + Dc,n = (20) 

with Dc — —0.0237 and D-ji — —0.0266. \Dc,r\ represents the distance separating the center 
of the circle {Y.c,n) to the origin of the poincare sphere. This proves that the planes are 
not going through the center of the sphere. It was checked after lengthy calculations that if 
\B\ = \C\ in the Jones matrix (see equation (1)) then D = 0. This shows that the property 
7^ is a characteristic of planar chirality (i.e, the condition \B\ ^ |C|). The radius of 
each circle {J2c,n) is given by rc,Ti — ^^(1 — Dj^ j^) and we have rc — 0.9997 and rc — 0.9996 
which are slightly smaller than r = 1 in agreement with the fact that Pc,n are not going 
through the center of the sphere. 
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